Classical Diffusion of a quantum particle in a noisy environment 
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We study the spreading of a quantum-mechanical wavepacket in a one-dimensional tight-binding 
model with a noisy potential, and analyze the emergence of classical diffusion from the quantum 
dynamics due to decoherence. We consider a finite correlation time of the noisy environment, 
and treat the system by utilizing the separation of fast (dephasing) and slow (diffusion) processes. 
We show that classical difi'usive behavior emerges at long times, and we calculate analytically the 
dependence of the classical diffusion coefficient on the noise magnitude and correlation time. This 
method provides a general solution to this problem for arbitrary conditions of the noisy environment. 
The results are relevant to a large variety of physical systems, from electronic transport in solid state 
physics, to light transmission in optical devices, diffusion of excitons, and quantum computation. 



Introduction.- The dynamics of quantum particles, de- 
scribed by the Schrodinger equation, manifests rich and 
exciting physics. The dynamics of a quantum particle in 
homogenous or periodic media is ballistic - i.e., the mean- 
square displacement of a quantum particle grows linearly 
with time, in contrast to the diffusive propagation of a 
classical particle. In many realistic cases, such quantum 
particles are exposed to time-dependent noisy environ- 
ments, and decoherence effects become important. In 
this work we focus on the emergence of seemingly classi- 
cal diffusion from the quantum dynamics under such con- 
ditions. Specifically, we calculate the classical diffusion 
coefficient as a function of the parameters of the noise - 
the magnitude and the correlation time. Our analytical 
calculations are found to be in excellent agreement with 
numerical simulations. 

To formulate the transition from quantum to classical 
diffusion, we use a tight-binding model with a stochas- 
tic potential. The problem of the spreading of a quan- 
tum wave-packet in a dynamic, time-dependent poten- 
tial has received much attention in the last three decades 
[3, S, S, i, S, H 0, Hi- More recently this problem was 
discussed in the context of quantum random walks with 
decoherence and quantum computation [lo| . This sim- 
plified theoretical problem is of great relevance to various 
experimental systems, such as particle diffusion in molec- 
ular crystals [3], diffusion of excitons I 111, pho ton prop- 
agation in coupled waveguide lattices [1^, ll3|, and de- 
struction of Anderson localization by nonlinearity or di- 
mensional crossover effects [il[i5l,[ll[l3- 
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The problem 

of dynamics and dephasing due to coup ling to a thermal 
bath is also closely related [l8, 1^, 20|. Early works [ij 
showed that unlike the static disorder case, where Ander- 
son localization sets in [21} , diffusive behavior takes over 
at long enough times for a tight-binding model with time- 
dependent disorder. Madhukar and Post extended 
this to the case of on-site as well as off-diagonal disorder. 
Various works [1, @, [2§| dealt with the problem on a con- 
tinuous lattice, which can show superdiffusive behavior. 
Heinrichs Q showed the correspondence between discrete 
(tight-binding) and continuous scenarios. All of these 
works assumed, for theoretical simplicity, delta-function 



correlations in time (white noise). The corresponding ex- 
perimental assumption for electron diffusion in molecular 
crystals, for example, is a temperature higher than the 
Debye temperature . Clearly, one would like to extend 
the theoretical understanding to more realistic cases of 
finite correlation times of the noisy environment. A first 
step in that direction was undertaken by [2^, using a 
perturbative analysis in the correlation time r. In this 
work we extend this to arbitrarily large correlation times. 
We use a novel method, relying on the separation of fast 
and slow processes (dephasing versus diffusion), to derive 
accurate analytical solutions for this problem. 

Model and derivation.- We consider the Schrodinger 
equation for a one-dimensional tight-binding model, tak- 
ing h ~ I. The equations governing the process are: 



dA, 
i—r- 
dt 



r(A,+i + A,_i) + e(j,t)Aj, 



(1) 



where Aj is the amplitude at site j and T is the nearest- 
neighbor tunneling. The noise term S,{j,t) is assumed to 
be Gaussian and uncorrelated between different sites, and 
will be characterized by a time correlation function C: 



(2) 



It is clear from the definition that C{—t) = C{t). The 
typical time of the decay of C{t) is defined as the cor- 
relation time r, and the typical magnitude is defined as 
W, i.e., C(0) = W'^. We shall assume that T -^W, from 
which follows that the dynamics of the phases of Aj+i 
and Aj^i will be driven mainly by the second term of 
Eq. ID). In the lowest order approximation A^ can be 
written as: 



|A°|e- 



(3) 



with</.,(0=(A,(0) + /o^(j,t')di'. 

Due to the random noise driving it, the phase e""^^*-*^ 
will also be a random variable, with a correlation function 
characterized by a correlation time (dephasing time) ^ 
T. We shall later show that this correlation function will 
be related to the diffusion constant of the particle. 
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The problem of finding the correlation function, de- 
fined as C4,{t - t') = {^pr'43(t)gi.<t>o(t'))^ is equivalent to 
that encountered in the physics of nuclear magnetic res- 
onance (NMR), where spins lose phase coherence due to 
the randomly fluctuating magnetic field created by the 



other spins, 
tain: 



Plugging in the solution for (pj (t) , we ob- 



(4) 



where we have omitted the site index from ^ (j, t) since 
this calculation is essentially for an independent, arbi- 
trary site. 

Formally, calculating this average is similar to the de- 
termination of Debye- Waller factors and yields: 



_lim 



This is an exact result for the correlation function of 
the dephasing process. A similar result was obtained in 
[ist - in the context of Dicke narrowing. When At ^ r, 
we can approximate the integral in the exponent of Eq. 
© by At C{t) = where /? is non-universal 

and depends on the form of C{t). For example, an ex- 
ponential decay gives /3 = 1 , while a correlation function 
decaying linearly to zero gives (3 = 1/2. Therefore in this 
regime: 



C^{At) « e 



(6) 



Thus for times large compared to r, the correlation 
time of the noise, the correlations of e~*'^ decay expo- 
nentially with a dephasing time ~ 1/W^t. This is 
reminiscent of the well-known NMR phenomenon of mo- 
tional narrowing (26| . 

For short times. At ^ r, we can approximate C{t) w 
C(0) = in the integral of Eq. and obtain: 



C4At) « e- 



(7) 



Remarkably, in this limit the dephasing is insensitive to 
the form of the correlation function, and is given by a 
Gaussian with T0 = \/2/W. 

Let us now proceed to analyze the dynamics of the 
probability distribution of the particle, characterized by 
the set of probabilities pj = l^jP- It is clear that with- 
out the tunneling term T the probabilities will remain 
constant, since the sites are uncoupled. However, due to 
the noise term, the phase of the amplitudes become a 
random variable, which, as we shall now show, will lead 
to classical diffusion of the probabilities. 

From Eq. ([T]), one obtains the exact relation: 



dt 



2TIm[A*Aj+i+A*Aj_i] . 



(8) 



Let us take an ensemble average of Eq. ([S]) . To the ze- 
roth order approximation, given by Eq. Q, the phase of 
1 are independent. This will lead, upon 



Aj,Aj^i and Aj 



taking the time-average, to the incorrect result (^) = 0. 
To proceed to the next order A^ in T/W, we have to solve 
the following first-order, linear differential equation: 



Defining an integration factor fij = 
this can be rewritten in the form: 



+ (9) 



d[A]^i,] 
dt 



T 



(10) 



Upon integration we obtain: 



T 



Jo 



fi,it')iA^^,+A^_,)dt'. (11) 



A similar equation holds for Aj+i and Aj-i. 
Plugging this into Eq. ^ and taking the ensemble- 
average we obtain: 



,dP^ 
' dt 



2rim 



(AO + Af)*(AO +Ai ) 



(12) 



where <-> denotes identical terms upon the substitution 

J + 1- > i - 1- 

A'j,A'j_^_i and j4°_j are independent, thus, only three 
terms remain in the average. The first of these is: 



/= (2T2Rc[-(A°) 



A^j+i(tO 
Mj+i(t) 



Aj"dt']), (13) 



while the other two are of similar form, obtained upon 
changing + and j ^ j ± 1. 

Let us evaluate /. Plugging in A^ = Pj(i)e-*(*), and 
the explicit form of /ij_|-i(t), we get: 



Jo 



gj[0,(*)-0,(i')]g»[03 + l(t)-</'j + l(t')]\ 



(14) 



Now comes the crux of the matter: the assumption 
T <C W, already utilized to allow a perturbative treat- 
ment in the leading order of T, also leads to a separation 
of timescales between the rate of change of the probabili- 
ties and the dephasing. This will be checked explicitly at 
the end of the calculation. Under this assumption, Pj{t) 
can be taken out of the averaging. Thus / = Pj{t)Q, 
with Q = (J* e'['^^'(*)"'^^'(*')le'['^J+i(*'"'^5+i(*')l). 

Assuming t 3> t^, we obtain that Q is a constant, 
related to the dephasing correlation function C0(AT) 
found before, through the relation Q = /g C^it)- 
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Using this result Eq. takes the form: 



dt 



2T^Q[P,+i + P,_i - 2Pj 



(15) 



where we have omitted the averaging notation. This is 
the well-known diffusion master equation. This implies 
the (ensemble-averaged) probabilities will indeed diffuse, 
with a diffusion coefficient D ~ 2T^Q. 

Combining this with Eq. ([5]) for the dephasing corre- 
lation function, we obtain our main result: 



D = 2T^ Cl{t')dt' = 2T^ / [f ~t)c{t)dt^^, ^ 

Jo Jo 

(16) 

For the case of short correlation time tW ^ 1, we can 
take C^{t') from Eq. and obtain: 



D 



rp2 



2T' 



(17) 



This complies with the results of Ovchinnikov [l|, who 
worked in the limit of r — > (delta-function correlations 
of the noise). 

The main advantage of our method, other than the 
clear physical picture gained for the diffusion process, 
is that we obtain explicit analytical results also for the 
regime of long correlation times, as well as the crossover 
between the two regimes. For long correlation times 
tW > 1, Eq. O gives C^{t'). Plugging this into Eq. 
(fTB)) we obtain: 



D 



w 



(18) 



The crossover between the two regimes (short and long 
noise correlation times), is also described by Eq. ([TB)) . 
The results depend on the form of the noise correlation 
fmiction C{t). Notice that for tW ^ 1 Eqs. dTT]) and 
(fT8|l give a result of the order of T^r, which is where the 
two limits 'match'. 

We have tested the theory numerically, on a finite lat- 
tice. The procedure is as follows: 

1. We solve Eq. ([1]) numerically, starting with a prob- 
ability distribution concentrated on a single site, and 
with a fluctuating disorder. For convenience, we take 
noise with linearly decaying correlations, i.e., C(t) = 
W^{1 - \t\/T)9{T - \t\), where 9 is the Heaviside step 
function. The wavepacket spreads, as is illustrated in 
Fig.IIl 

2. We ensemble- average the probability distribution 
over many realizations of the noisy environment. We 
find that the standard deviation of the ensemble-averaged 
probability distribution grows as a{t) ~ \/2Dt, as pre- 
dicted by the theory. This is demonstrated in Fig. [51 



From this dependence we extract the diffusion constant 
D. We find indeed that D depends on the strength of 
the noise W as well as its correlation time r. 

Fig. [3] compares the result of the simulations with the 
theoretical prediction of Eq. p6|) . with a linearly decay- 
ing C(t), as used in the simulations. Since the diffusion 



constant is proportional to T 



know from dimen- 



sional analysis that D = T Tf{WT), where f{x) is de- 
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FIG. 1; The spread of a quantum-mechanical wavepacket in 
a perfectly periodic lattice (T = 1, = 0) showing ballistic 
spreading (a,b), versus a noisy lattice with T = 1, t — 0.01, 
W — 20 showing diffusive spread (c,d). The wavefunction 
is confined to a single lattice site at t=0. The probability 
distribution of the quantum particle after some propagation 
is plotted on a semilog scale in (a) and (c) , showing a ballistic 
profile and a diffusive (gaussian) profile, correspondingly. 




Simulation results 
Theory with D=2T^/W^t 



2000 4000 6000 8000 10000 

Time 

FIG. 2: A plot of the width versus time of a quantum me- 
chanical wavefunction spreading in a noisy potential , showing 
difi'usive propagation. The simulation parameters are T = 1, 
r — 0.01, W — 20. The data is averaged over 100 realizations 
of the disordered potential. The numerical results (open cir- 
cles) fit the theory (line), with no fitting parameters. 
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tcrmincd from Eq. p6|l . This allows for a unified presen- 
tation of both regimes, which are manifested as different 
asymptotic regimes of the function f{x): for Wt <^ 1 
we know from Eq. pT]) that f{x) ~ l/x^ with a non- 
universal proportionality constant, while for Wt 3> 1 we 
have f{x) = ^/2tt/x, as deduced from Eq. (fT^ . Indeed, 
we numerically obtain scaling of this form, where the 
data for different runs collapse onto a single curve, when 
scaled correctly. Fig. [3] shows f{x) derived numerically, 
and compares it to the prediction of Eq. (fT6|) . Note that 
there are no fitting parameters. The excellent fit confirms 
the validity of the analytical approach used. 

Summary.- We have considered a tight-binding model 
with random on-site energies, fluctuating in time with 
an arbitrary correlation time. Using the separation of 
timescales between the (slow) diffusion process and the 
(fast) dephasing, we managed to map the complex evo- 
lution of the probabilities onto the problem of dephasing 
of a single site. For a noisy environment with short cor- 
relation time (compared to the reciprocal strength of the 
disorder, l/W), the phenomenon of motional narrowing 
is manifested in the diffusion constant. In the regime of 
long correlation time, we show that the diffusion constant 
is inversely proportional to the strength of the disorder. 
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FIG. 3: A fit of tfie theoretical results to numerical calcu- 
lations. The diffusion constant was found numerically, for 
correlation times r = 0.01 (triangles), r — 0.1 (squares) and 
T = 1 (circles) and noise magnitudes ranging from W = 2 to 
W = 20. Each point in the graph is an average over 50 real- 
izations. As expected theoretically, there is a crossover from 
a regime with slope —2 (i.e., dependence), to a regime 

with slope —1 ( 1/x dependence). The different curves over- 
lap at points with the same value of Wt. The form of the 
crossover is given by Eq. (|16p . integrated numerically for the 
intermediate regimes. No fitting parameters are used in the 
comparison. 



We also give a solution for the crossover regime. Our re- 
sults show excellent agreement with simulations, inadver- 
tently confirming previous results obtained for the some 
particular cases studied in the past, but continuing much 
beyond to give a general solution for the full range of the 
parameter space. This includes quantum systems under 
more realistic noise environments conditions, which have 
not been studied before. Our results and exact solutions 
are applicable to various quantum systems, from diffu- 
sion of electrons and excitons to photon propagation in 
coupled waveguides, and to the general question of quan- 
tum random walks with decoherence. 
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for useful discussions. A. A. acknowledges funding by the 
Israel Ministry of Science and Technology via the Eshkol 
scholarship program. Y.L. acknowledges support by the 
Israel Academy of Science and Humanities via the Adams 
scholarship program. 
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